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Abstract. We present an analysis of low-resolution infrared spectra for 
25 brown dwarf candidates in the NGC 1333 molecular cloud. Candidates 
were chosen on the basis of their association with the high column density 
cloud core, and near-infrared fluxes and colors. We compare the depths 
of water vapor absorption bands in our candidate objects with a grid of 
dwarf, subgiant, and giant standards to determine spectral types which 
are independent of reddening. These data are used to derive effective 
temperatures and bolometric luminosities which, when combined with 
theoretical tracks and isochrones for pre- main sequence objects, enable 
us to estimate masses and ages. Depending on the models considered, a 
total of 9 to 20 brown dwarfs are identified with a median of age of <1 
Myr. 



1. Introduction 

The NGC 1333 reflection nebula and its associated dark cloud L 1450 are part 
of a chain of molecular clouds in the Perseus region (Sargent 1979; Loren 1976). 
Analysis of Hipparcos data suggests a distance of 300 pc (de Zeeuw, Hoogerwerf, 
& de Bruijne 1999; Belikov et al. 2002). The observations of emission-line stars 
and Herbig-Haro objects first established NGC 1333 as an active region of star 
formation (e.g., Herbig 1974). Surveys of the cloud at near-infrared wavelengths 
have revealed a large population of low mass stars (Strom, Vrba, &i Strom 1976; 
Aspin, Sandell, & Russell 1994; Lada, Alves, & Lada 1996; Wilking et al. 2002). 
The young stellar objects are concentrated into a northern and southern cluster, 
each with about 70 members. Identification of the lowest mass objects in the 
NGC 1333 cloud, as well as estimates for their ages and masses, requires infrared 
spectroscopy. Ultimately, the population of very low mass objects will be used 
to define the shape of the mass function and determine the number and total 
mass of brown dwarfs in the young cluster. 
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2. Source Selection and Observations 

Brown dwarf candidates were selected by comparing their positions in a near- 
infrared color-magnitude diagram with theoretical models for pre-main sequence 
objects by D'Antona & Mazzitelli (1997, 1998, hereafter DM97). An age of 1 
Myr was assumed (Lada et al. 1996). Near- infrared surveys of the northern 
cluster (Wilking et al. 2002) and the southern cluster (Aspin et al. 1994) were 
used to construct the K vs. (H-K) diagram. The majority of the 25 candidates 
lie in a region of the diagram where M<0.1 Mq, A„ <10 mag, and K<14 mag. 

Infrared spectroscopic observations were accomplished using the 3-m NASA 
Infrared Telescope Facility at Mauna Kea, Hawaii in 2000 November 10-13. The 
observations were aided by an Internet2 link between UM-St. Louis and Mauna 
Kea. Spectra in the K band were obtained for 25 brown dwarf candidates, eight 
M giants, and eight M subgiants using the 256 x 256 InSb facility infrared camera 
(NSFCAM) with the HKL grism and a 0.3" pixel" 1 scale. The 0.6" slit provided 
a resolution of R = A/AA ~300 over the 2.0 ^m-2.5 /Ltm band. 

3. Results 

3.1. Spectral Classification Using a Water Vapor Index 

In the K-band, broad water vapor absorption bands centered at 1.9 /im and 2.5 
fim provide the best means for deriving spectral types for faint M stars. These 
bands are not only sensitive to temperature for cool stars (Jones et al. 1994) but 
also are well-resolved in low resolution spectra. Water vapor indices have been 
used to derive spectral types for very low mass young stars in the p Ophiuchi 
cloud (Wilking, Greene, & Meyer 1999, hereafter WGM99; Cushing, Tokunaga, 
& Kobayashi 2000), IC 348 (Najita, Teide, k Carr 2000), the Orion Nebula 
cluster (Lucas et al. 2001), and the Taurus cloud (Itoh, Tamura, & Tokunaga 
2002). In WGM99, we defined a Q index that was independent of reddening as 

Q = (F1/F2)(F3/F2) 1 - 22 (1) 

where Fl, F2, and F3 are the average fluxes in narrow bands covering 2.07-2.13 
/im, 2.2695-2.2875 pm, and 2.40-2.50 ^m, respectively. 

A linear fit to a plot of the Q index vs. the optically-determined M spectral 
type for dwarf standards is shown in Fig. 1 for stars from M0.5V-M9V and 
yields the following relation 

MV subclass = (-20.15 ± 0.48) x Q + (18.49 ± 0.55) (2) 

with a correlation coefficient of r=0.98 for a sample of 12 measurements. 

3.2. The Water Index as a Function of Surface Gravity 

A sample of eight M subgiants were selected for infrared spectroscopic observa- 
tions from the MBM 12 and IC 348 young stellar populations. The median ages 
of these populations are estimated to be 2 Myr. As shown in Fig. 1, the linear 
fit to the subgiant sample is consistent with the dwarf relation and there is no 
apparent sensitivity of the Q index to surface gravities in the dwarf to subgiant 
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Water Index vs. M Spectral Type 




Figure 1. The Q index vs. M spectral type for dwarf standards 
(squares), subgiants (diamonds), and giants (triangles). Linear fits to 
the dwarf and subgiant data are nearly identical. 



range. This is in agreement with the findings of Jones et al. (1995) who, in 
comparing infrared synthetic spectra, saw little dependence on the depths of 
the water lines with surface gravity between log(g)=4.0 to 5.0. The Q values 
for the 8 giants observed show a much greater scatter from a linear fit. In our 
sample, dwarfs always display deeper water vapor absorption than giants for a 
given spectral type. However, the Q index for giants can be smaller than that 
of dwarfs due to enhanced absorption by CO in the 2.40-2.50 /im band. Hence, 
the Q index is not a useful measure of water vapor absorption for objects with 
the surface gravities of giant stars. 



3.3. Spectral Types 

Twenty-four brown dwarf candidates displayed late-type photospheres with ab- 
sorption due to water vapor and CO characteristic of high surface gravity. Since 
the surface gravities of objects in this sample are likely to resemble those of 
subgiants, spectral types were derived using Eqn. 2 which is appropriate for 
both subgiants and dwarfs. Spectral classifications using the water vapor index 
relation reveals spectral types between M2.4 and M8.3 with typical uncertainties 
of ±0.5 subclass. Seventeen of 24 candidates displayed spectral types >M6.0. 
The visual extinctions implied by these spectral types are typically 2-3 mag 
and all less than 10 mag. The majority of the candidates appear to be young 
objects of very low mass. The projection of this sample on the high column 
density core (A„ >10 mag) coupled with the relatively low visual extinctions 
of the candidates minimizes the chance that some objects are background M 
stars. The detection of x-ray emission from 15 of the 24 candidates by Getman 
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Figure 2. H-R diagrams for the brown dwarf candidates using the 
DM97 models (left panel) and the Burrows et al. models (right panel). 
'X's mark the positions of x-ray emitting candidates. 



et al. (2002), including 3 tentative detections, argues in favor of the youth of 
this sample, as x-ray luminosity and stellar activity are known to decline with 
age for K and M stars (Fleming, Schmitt, & Giampapa 1995). 

3.4. Masses and Ages of the Brown Dwarf Candidates 

We followed the procedure outlined in WGM99 to estimate the masses and 
ages of the brown dwarf candidates. Effective temperatures computed for each 
brown dwarf candidate have typical uncertainties of ±95 K. However, systematic 
effects could affect the derived temperatures. The presence of moderate veiling 
in 3 objects could lead us to overestimate their effective temperatures by ~150 
K. The assumption of dwarf, rather than subgiant, surface gravities could also 
systematically affect our results. For spectral types of M2 and later, giant star 
temperatures are warmer than dwarfs by 300-500 K for stars of the same spectral 
type (e.g., Fig. 4 in Itoh et al. 2002). Hence by assuming a dwarf surface gravity, 
we may be underestimating the effective temperature of an M2 subgiant by ~150 
K and of an M6-M9 subgiant by ~250 K. 

In Fig. 2, we have plotted our brown dwarf candidates on Hertzsprung- 
Russell diagrams overlaid with the theoretical tracks and isochrones from the 
models of of DM97 and Burrows et al. (1997). An H-R diagram using the Baraffe 
et al. (1998) models is presented in Wilking et al. (2002). For the DM97 models, 
mass estimates range from <0.02 Mq to 0.25 Mq, with 16 objects at or below 
the hydrogen-burning limit. The median age for the sample is 0.3 Myr. If a 
temperature scale intermediate between dwarfs and giants is used, then these 
objects should be shifted to higher temperatures by 250 K and only 9 would 
remain as brown dwarfs. The Burrows et al. models yield lower masses and 
younger ages with 23 objects at or below the hydrogen-burning limit and 20 of 
these remaining as brown dwarfs if shifted to higher temperatures by 250 K. The 
median age for the sample implied by the Burrows et al. models is 0.1 Myr. 
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3.5. Surface Gravity Estimates 

Surface gravities (g=GM/R 2 ) can be estimated for the brown dwarf candidates 
using the DM97 mass estimates and radii derived from the luminosity and ef- 
fective temperature. Values range from log(g)=3.0 to 4.5 in cgs units with a 
median value of 3.3. These estimates assume a temperature scale, intrinsic col- 
ors, and bolometric corrections derived from dwarf standards and only increase 
under the assumption of giant properties. Therefore, the surface gravities of 
our objects more closely resemble dwarfs than giants, implying that our use of 
a dwarf temperature scale has led to an underestimation of the temperature for 
M6-M9 stars of <250 K. When taken along with the results of the two sets of 
models, we conclude that at least 9 of the candidates can be identified as true 
brown dwarfs. 
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